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Abstract 


Results of continuing efforts toward validating the buildup-factor method and the baryon 
transport code (BRYNTRN), which use the deterministic approach to solving radiation 
transport problems and are the candidate engineering tools in space radiation shielding analyses, 
are presented. A simplified theory of proton-buildup factors assuming no neutron coupling 
has been derived to verify a previously chosen form for parameterizing the dose-conversion 
factor that includes the secondary-particle buildup effect. Estimates of dose in tissue made by 
the two deterministic approaches and the Monte Carlo method are compared for cases with 
various thicknesses of shields and various types of proton spectra. The results are in reasonable 
agreement but there is some overestimation by the buildup-factor method when the effect of 
neutron production in the shield is significant. Future improvement, including neutron coupling 
in the buildup-factor theory, is suggested to alleviate this shortcoming. Impressive agreement 
for individual components of doses, such as those from the secondaries and heavy-particle 
recoils, is obtained between BRYNTRN and Monte Carlo results. 


Introduction 

As NASA continues to develop a vigorous space program, tools for the analysis of optimum 
shielding against space radiation are a continuing requirement. The tools must ultimately 
account for the complex mixtures of radiations in the space environment and be complete in 
the description of the physical processes involved to minimize engineering-design error. Still, a 
complete model must be computationally efficient and numerically accurate to be a useful tool 
for design work. For the model to be used with confidence, efforts toward model validation 
must be made. 

Monte Carlo computer codes have been written that meet many of the above requirements 
(ref. 1). However, the enormous computational time and storage requirements have impeded 
their usefulness in the space program. Other alternatives, such as BRYNTRN, a nucleon 
(baryon) transport code (ref. 2), and a galactic cosmic ray (GCR) transport code (refs. 3 and 4), 
both developed at Langley, use the deterministic approach, which imposes fewer computational 
requirements. In reference 2, comparisons were made between BRYNTRN and Monte Carlo 
calculated doses in bare (no shielding) tissue exposed to incident monoenergetic protons, and 
the results were in reasonable agreement. 

Current efforts at Langley include (but are not limited to) pursuing further improvement 
and validation of the transport codes and developing a simple, easy-to-use, application-oriented 
code for estimating doses that result from incident space protons (refs. 5 and 6). The latter 
effort makes use of the buildup factors, which compress the detailed computation of secondary- 
particle buildup into the dose-conversion factor. The buildup-factor method is clearly the most 
computationally efficient method if one is only interested in doses as the end result, but it is 
limited to materials for which buildup factors are known. 

The purpose of this study is to continue the comparison effort made in reference 2, to 
include dose estimation in tissue behind shielding materials, and to extend the comparison to 
the newly developed buildup-factor code, BRYNTRN, and the Monte Carlo method. Dose 
calculations are made for protons with either continuous or discrete spectra where low-energy 
Monte Carlo or experimental data are available. Doses calculated for the penetrating solar flare 
of the February 1956 event by the buildup-factor method and BRYNTRN are also compared 
to ensure that the parameters used in the buildup factors are valid for the high-energy region. 
A detailed discussion of a simplified buildup-factor theory is given. 



Symbols 


a., 

B 

B A 

b 

D(x, E) 

Di 

D s (t s , x,E) 

E 

E r 

E, 

f(E,E') 

fir , r') 

HZE 

m 

P 

Q 

r(E) 

S(E) 

Si 
t . , 


X, Z 
0 

f(r) 

a 


4>{x, E) 
^(x,r) 
Subscripts: 
0 

P 

s 


parametric constants 

buildup factor or ratio of total dose to primary proton dose 
build up- factor correction 

slope parameter in secondary proton production cross section, cm 4 /g 2 

total dose (or dose equivalent) at depth x in exposed medium to inci- 
dent proton of energy E , rad (or rcm) 

ith order of dose (or dose equivalent) contribution, rad (or rem) 

total dose (or dose equivalent) at depth x in exposed medium behind a 
shield of thickness t s to incident proton of energy E , rad (or rem) 

proton energy, MeV 

= e (ro — x), residual proton energy after penetrating a distance x, MeV 
= e (r 0 - t s ), MeV 

production cross section for secondary proton of energy E by primary 
proton of energy E\ cm 2 /g-MeV 

= S(E)S(E,E f ), cm 4 /g 2 

high-energy heavy ion 

average number of secondary protons produced per nuclear event 
exp(— r) 
quality factor 

range of proton with energy E , g/cm 2 

stopping power of proton with energy E, MeV-cm 2 /g 

2 th order spectral average of stopping power, MeV-g/cm 2 

equivalent thickness of shield relative to tissue such that stopping 
powers of shield and tissue are equal, g/cm 2 

distance or range, g/cm 2 

parameter in secondary proton production cross section, cm 2 /g 
energy of proton with a range r, MeV 
macroscopic interaction cross section, cm 2 /g 
optical thickness for proton 

proton fluence spectrum at distance x, protons/cm 2 -MeV 

= S(EWx,E), g- 1 

initial value 
proton 

shield material 
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where m is the average number of protons produced per nuclear event. Although m and a are 
in reality functions of E f , our current interest is m monoenergetic boundary conditions as 

if>(Q,E) = 6(E-Eo) (8) 


When m and o are evaluated at the beam energy Eq, the corresponding boundary condition 
on 0 is 

“0(0, r) = 6{r - r 0 ) (9) 

The high-energy production cross section is an exponential function of E f — E and is used to 
approximate equation (4) as 

f(r, r) = bcxp[—/3(r f — r)] (10) 

The normalization in equation (7) requires 


b = 


ma0 

1 - exp (-/?r 0 ) 


( 11 ) 


and 0 = 0.01 cm 2 /g. Equation (5) may be solved by perturbation theory (ref. 2) to obtain 


oo 

tp{x,r) = 

2=0 


( 12 ) 


where 



ipo(x, r) = exp(— crx)6(r + x — tq) 

(13) 


f x ~ 

ipl(x. r) = exp(— ax) / f(r + z, tq — x + z)dz 

Jo 

(14) 


rx roc 

0 n+ l(x,r)= / dzexp(-az) / dr'f(r + z, r f )'ip n (x — z, r f ) 
Jo Jr+z 

(15) 

Equation 

(12) may be reduced using equation (10). For example, 



01 (^i r) = x exp(—ax)f(r, tq — x) 

(16) 


1 2 

02 (x, r ) — 2 X oxp(—ax)b • (ro — x — r) • /(r, tq — x) 

(17) 

The successive contributions to dose may now be calculated as follows: 



roc 

Dq(x) = / exp(— ax)S(r + x — ro)dE = 5[c(ro — x)] exp(— ax) 
Jo 

(18) 


D\{x) = 5i[e(ro — x)]max exp(-ax) 

(19) 


- o 2 x 2 

D 2 (x) = 52[e(ro x)]m 2 exp( ax) 

(20) 

where S\ and 52 are spectral averages of stopping power in which 5i(e) = 
5 2 (c) = 0(e 3 ) for small values of e. The total dose is then 

0{ e 2 ) and 


00 i - 

D(x) — S[e(ro - x)] exp(-cjx) + -{max) 1 exp(-(rx)S'j[f(ro — x)] 

2=1 *' 

(21) 
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Simplified Theory of Proton-Buildup Factors 

In passing through material or tissue, energetic protons interact mostly through ionization 
of atomic constituents by the transfer of small amounts of momentum to orbital electrons. 
Although the nuclear reactions are far less numerous, their effects are magnified because of 
the large momentum transferred to the nuclear particles and the struck nucleus itself. Many 
of the secondary particles of nuclear reactions are sufficiently energetic to promote similar 
nuclear reactions and to cause a buildup of secondary radiations. The description of such 
processes requires solution of either the integral or integrodifferential form of the Boltzmann 
equation as is done, for example, in the Monte Carlo method or BRYNTRN. To bypass such 
extensive calculation for the physical processes, the effect of secondary- particle buildup was 
folded into dose-conversion factors that relate the primary monoenergetic proton fluence to 
dose or dose equivalent as a function of distance in a slab of tissue or material (refs. 5, 7, 
and 8). 

A parameterization of the conversion factors, including the secondary-particle buildup in a 
slab of tissue, was first introduced by Wilson and Khandelwal (ref. 7). Their method allows 
reliable interpolation and extrapolation of known values obtained from Monte Carlo results. 
A refinement to the parameterized form of Wilson and Khandelwal was subsequently made 
(ref. 8). The work has also been extended to the tissue behind an aluminum shield, and an 
aluminum-buildup factor relative to tissue has been introduced (ref. 8). To elucidate and verify 
how these parameterized forms were obtained, a simplified theory that assumes no neutron 
coupling from aluminum is presented in the following sections and serves as a basis for future 
development. 

Buildup in Exposed Medium 

The Boltzmann equation for proton transport in the straight ahead approximation is given 

roc 

<j>{x,E)= f{E,E')(t>{x,E')dE' (1) 

Je 

where S(E) is the proton stopping power, a is the macroscopic interaction cross section, which 
we presently take as energy independent, and f{E,E f ) is the production cross section for 
secondary protons (ref. 2). Using the definitions 



rE 

r= / dE'/S{E') 

Jo 

(2) 


= S(E)cf)(x, E) 

(3) 

and 

f{rJ) = S{E)f{E,Ef) 

(4) 

allows equation (1) to be rewritten (ref. 2) as 


^(x,r)=e ax ip(0,x + r) 

rx roc 

1 + / dzexp(-az) / dr ! f{r + z y r / ) j 0(a: — z, r) 

Jo Jr+z 

(5) 

where the boundary condition is 

^(0, r) = S(E)<f>(0,E) 

(6) 


The secondary-particle production cross section is normalized as 

f E ' 

/ f(E, E f )dE = m<r (7) 

Jo 


f x -6E S ^ + ° 
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The dose buildup factor defined as the ratio of the t« ,tal dose to the primary proton dose is 
then 


OO , , 

E h(m rx -'SMro-x)] 


B(x , Eq) = 1 + 


i=l 


5[e(r 0 - x)] 


with the property that 


lim B(x, Eq) — 1 

x->ro 


which follows here from the neglect of the coupling between the proton and neutron fields. 
Wilson and Khandelwal (ref. 8) assumed that the buildup factor had the form 

B(x, E 0 ) = (ai + a 2 x + a 3 x 2 ) exp(-a 4 x) 


(22) 

(23) 

(24) 


where a 4 was chosen to satisfy equation (23). The dioice of a 4 is not to be governed by the 
nuclear cross section, but rather by the result that S t (e) ~ 0(e l+1 ) for small values of e. Tins 
conclusion is modified because of the presence of neutron production in the medium. 

The buildup-factor parameters are shown in table 1 and are to be used in the following 

l0mK D(x,Eq) = B(x,E 0 )Q[S(E r )}P(Eo)S(E r )/P(E r ) (25) 

where E r is the residual energy at r(E 0 ) - x.Q(S) is the quality factor as a function of the 

stopping power, and . 

P(E) = exp[-r(E)] ( 26 ) 


with 


r(E) = f E a(E')dE'/S(E') 

J o 


(27) 


Equation (27) expresses the full energy dependence of the nuclear cross sections. The optical 
thickness t(E) is given in table 2. Note that Q(S ) is set to unity when expressing absorbed 

dose. 

Buildup in an External Shield 

The effect of buildup in an external shield on the exposed medium (tissue) doses is now 
considered. If it is assumed that an equivalent distance in the shield can be defined such that 
the stopping power in equivalent distance units of the shield and exposed medium (tissue) are 
equal, then the Boltzmann equations of the two media differ only in their nuclear cross sections 

9 9 — <Kx,E)= / fs{E,E')4>{x,E')dE' 

J E 


dx dE 

d_ _ d_ 
dx dE 


S{E) + <7 S 
S(E) + a 


IE 

roc 


roc 

4>(x,E)= / f{E,E')4>{x,E')dE' 

J E 


(28) 

(29) 


For a monoenergetic beam on the boundary of the shield, the solution is given as 

i/>o(x,r) = exp(-cr.,x)<5(ro - x - r) 
rp\{x,r) = xexp(-<T s x)/.,(r,ro - x) 


(30) 

(31) 


The particles appearing at the media interface provide the boundary condition of the exposed 
medium; thus, 

V> 0 (0, r) = exp(— <7 s t s )£(r 0 -t s -r) + t s exp(-cr s t s )f s (r, H) -*») + ■• ■ ( 32 ) 
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where t s is the equivalent thickness of .he shidd. To evaluate the proton field in the exposed 
media, we may use equation (5) and the boundary value in equation (32) to obtain 

W{x, r) = exp(— < tx) exp(-a.,t.,)6(r 0 - t' x - r) + exp(-ax)<. s exp(-a,t s )f s {r + x, r 0 - f.,) 

+ exp(-a s f. s )xexp(-crx)7(r,ro - t s - x) + . . . (33) 

Utilizing equation (10), equation (33) may be rewritten as 

t/>(x, r) = cxp( tXgttf - ax) 6 (ro - t s - x - r) 

+ [tsfs(r, ro - t s - x) + xf(r, tq - t s - x)] exp(-a s t s - ax) + . . . (34) 

Similar to equation (21), we can write the total dose at depth x in an exposed medium behind 
a shield of equivalent thickness t s as 

D s (t s , x, Eq) = exp (—a s t s — rrx)S'[e(ro — t s — x)] T {m s cr s t s S l s i [e(ro — t s — x)] 

+ maxS 1 [e(r 0 - t„ - x)] } exp(-a x t s - ax) + . . . (35) 

Equation (35) may be written as 

D s (t s ,x, Eo) = D(t s -f x , Eo) 4- [exp(— a s t s ) - exp(-cr£ s )] D(x , E s ) 

+ f s exp(-<7x) {rrisCTsSsi [e{r {) - t s - x)} exp (~a s t s ) 

- rnaS[ [e(r 0 - t s - x)] exp {~at s ) } + . - - 

~ E>(t s + x, Eq) -f t s {(a - cr s )D(x, E s ) + exp(-crx)(m 5 a s - me r) 
x Si[e(r 0 ~t s -x)}} + ... ( 36 ) 

where E s = e(ro - t s ). The shield buildup factor can be defined relative to the exposed medium 
as 


X y Eo) — &s(tsf *£? Eo)/ D(t s + X, Eq) 

= 1 + ts[{<r - <T S )D(x, E s ) + exp (-ax)(m s a s - ma)Si]/D(t 9 + x, E 0 ) + . . . (37) 

Clearly, the coefficient of t s in equation (37) is reduced if <y — & s and reduces to a smaller 
contribution as m — ► This reduction results because the spectral distribution functions 

have 0 S = 0 for all materials, where 0 is largely determined from the proton-proton scattering 
amplitude. 

For the present work, the buildup- factor forms from reference 5, which possess the same 
quality as equation (37), are used and are given as 

D u , 0.02 Et s 

Eo) = 1 -f 1 E exp(-Q.Q22 t s ) (38) 

for the dose equivalent and 

B A (t s ,x,E 0 ) = 1 + 6(1+ ^ ) exp(-0.0U s ) (39) 
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for the absorbed dose, where E is in MeV and t s in gm/cm 2 . Equations (38) and (39) are used 
in the expression 

D s (ts, *, E 0 ) = B A (t s , x, Eo)B{t a + x, E 0 )Q[S(E r )}P(E)S(E r )/P(Er) (40) 

for the dose in the tissue behind an aluminum shield where Q is set to unity for the absorbed 
dose. Equations (38) and (39) are in agreement qualitatively with the derived form (eq. ( ) 

in that values of B A {t s ,x,E) are linear in t, for small values of t s and linear m E for sma 
values of E. Also, these equations are chosen to fit existing data that exhibit plateaus for large 
values of E or t s . Although the derived form (eq. (37)) has resulted from the ^sumption ° 
no neutron coupling, the chosen form with the constants as given in equations (38) and (39) is 
obtained based on data that include neutron coupling. 


Results and Discussion 

Since the buildup factors are functions of both energy and thickness, the first step in verifying 
such a method is to compare the results at various fixed (discrete) energies. In reference 2, 
comparisons were made between BRYNTRN and Monte Carlo results for monoenergetic protons 
at various energies, and they were in reasonable agreement considering the numerical difficulty 
involved in discrete energy calculations (ref. 2) with BRYNTRN. For the buildup-factor method, 
comparisons made between Monte Carlo and experimental data are shown m figures 1 to 3. 

The dose and dose equivalent calculated as functions of depth in tissue with and without an 
aluminum shield are shown in figures 1 and 2 for normal incident protons at discrete energies 
of 400 660 730 1500. and 3000 MeV. The limited Monte Carlo results with 0 g/cm shielding 
are obtained from reference 9. The calculated values with the buildup-factor method are in 
reasonable agreement despite the crudeness in the buildup parameters chosen. Although there 
are no Monte Carlo data available with shielding at discrete energy, the doses calculated with 
30 g/cm 2 of aluminum shield by the buildup-factor method are also presented in the figures or 
qualitative comparison. In general, the dose is increased because of the presence of the shield. 
The increase in dose over those with no shielding results from neutrons producec in a uminuni, 
especially in the first few centimeters of the tissue. For protons at the lowest energy (400 MeV), 
the Bragg peak appears at 55 cm in depth as the protons approach their limiting range while 

traveling through the shield and tissue. c 

Figure 3 shows the comparison of the experimental data (fig. 4 of ref. 7), the buildup- factor 
calculation, and the interpolated Monte Carlo results (refs. 9, 10, and 11) for the absorbed 
dose in tissue that is exposed to a proton beam of 592 MeV. Also shown are the earlier (re . ) 

buildup-factor calculations for uncollided primary and total absorbed dose. The dose from 
the buildup, which is the difference between the total and uncollided primary, is substantial. 
The usual Bragg peak is also obvious for both the analytical and experimental resu ts e 
buildup-factor calculations are in general within the uncertainties of interpolated Monte Carlo 
values and are in reasonable agreement with the experimental data. 

To verify both BRYNTRN and the buildup-factor method in case of a continuous energy 
spectrum of incident protons, dose calculations were made (figs. 4 and 5) for a shielded tissue 
exposed to a typical solar-flare spectrum for which Monte Carlo results were availa e (re s 
and 13). The flare spectrum taken from references 12 and 13 of the Webber form (ref. ) 
is exponential in rigidity with characteristic rigidity P 0 = 100 MV and is normalized to 
10 9 protons/cm 2 with energy greater than 30 MeV. Only the portion of the spectrum between 50 
and 400 MeV was considered for the Monte Carlo calculation (refs. 12 and 13). Nevertheless, or 
the current calculations, the high-energy cutoff at 400 MeV was ignored; very small differences 
of a few percent were found (figs. 4(a) and 5(a)), because the spectrum contains very few highly 
penetrating energetic protons, which may become significant only at depths beyond current 
interest. The aluminum shields (fig. 4) and the iron shields (fig. 5) for the tissue were 20 g/cm 
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in thickness. The buildup-factor results presented for iron, however, were obtained with the 
same values of buildup parameters chosen for aluminum, because the nuclear-reaction cross 
sections are roughly the same for both materials at the energies of interest (below 400 MeV). 
The total doses by the buildup- factor method are in good agreement with both Monte Carlo and 
BRYNTRN results (see figs. 4(a), 5(a), and 5(b)) except at the first few centimeters in tissue 
depth, where the effect of neutron production in the shield is significant. Therefore, future 
improvement is needed in the buildup-factor method for treating neutron coupling explicitly. 

Detailed dose comparisons are also given between BRYNTRN and Monte Carlo results in 
figures 4(b), 4(c), and 5(c) to 5(f) for the aforementioned low-energy spectrum. The primary 
doses are in good agreement where, for the worst case, the discrepancies came partly from the 
intrinsic nature of the Monte Carlo fluctuations. This fluctuating phenomenon is even more 
severe for the total secondary dose (figs. 4(c), 5(e), and 5 (f)). (The total secondary dose is the 
sum of the secondary- proton and secondary-neutron doses, as they are termed in references 12 
and 13.) The heavy-ion-recoils dose and dose equivalent of BRYNTRN show that the actual 
physical dose from heavy recoils may not be important, but that the dose equivalent can be 
significant because of the large quality factor. 

Dose calculations are also made for a continuous spectrum that contains more high-energy 
protons. Since there are no Monte Carlo results available in the high-energy range, the February 
1956 solar-flare event is chosen for comparison between BRYNTRN and the buildup-factor 
method. The solar-flare spectrum given as the integral fluence form in protons/cm 2 is 


<M> E) = 1.5 x 10 9 exp +3 * 10 8 exp (-^ 3 ^°° ) (41) 

where E is the energy in MeV. The results shown in figures 6 (a) and 6 (b) are for the dose and 
dose equivalent in tissue with 0, 10, and 30 g/cm 2 of aluminum shield. The agreement between 
these two deterministic methods is reasonably good. With the future improvement to include 
neutron coupling, the buildup-factor method, which is gaining in computational efficiency for 
flare-dose analysis, would probably be favored in parametric studies of spacecraft shield design. 

Concluding Remarks 


Comparisons between the buildup-factor method, the baryon transport code (BRYNTRN) 
and the Monte Carlo method are made for the calculated doses in tissue behind various 
thicknesses of shielding with exposure to various proton spectra. The results are in reasonable 
agreement, but there is some overestimation by the buildup factors when the effect of neutron 
production in the shield is significant. Future improvement, including neutron coupling in 
the buildup-factor theory, should alleviate the shortcoming. Impressive agreement for various 
components of doses is obtained between BRYNTRN and the Monte Carlo results. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
August 22, 1990 
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Table 1. Buildup-Factor Parameters for Tissue 


Dose equivalent 


GeV 

a \ 


^3 

<24 

ai 

a 2 


04 

).03 

1.00 

0 



0 

1.00 

0 



0 

*.06 

1.20 

0 



.0130 

1.07 

.010 



.0100 

.10 

1.40 

.020 



.0300 

1.10 

.040 



.0260 

*.15 

1.45 

.065 



.0385 

1.12 

.060 



.0310 

.20 

1.50 

.080 



.0400 

1.15 

.062 



.0320 

.30 

1.60 

.100 



.0330 

1.20 

.068 



.0260 

.40 

1.80 

.120 

J 


.0228 

1.24 

.071 

> 


.0228 

.73 

3.00 

.140 

.00031 

.0160 

1.40 

.090 

.00010 

.0150 

1.20 

3.90 

.150 

.00130 

.0150 

1.67 
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Interpolated values. 


Table 2. Total Tissue Optical Thickness for Protons 


E . GeV 


E . GeV 

.300 
.500 
.700 
.000 
.200 
.400 
.600 
.800 
.000 
.000 
.000 
.000 
.000 
.500 
.000 


t ( E ) 

6.57 

8.03 

9.52 

11.76 

13.27 

14.78 

16.29 

17.79 

19.29 
26.62 
33.81 

40.84 
47.75 
57.91 

67.85 
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Figure 3. Comparison of predicted dose in tissue with experimental data for normal incident 
protons at 592 MeV. 
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Figure 4. Doses in tissue behind 20 g/cm 2 of aluminum shield to normal incidence of a Webber 
solar-flare proton spectrum. 100-MV rigidity. 
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(b) Primary proton dose. 
Figure 4. Continued. 
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(c) Total secondary and heavy-ion recoil dose. 
Figure 4. Concluded. 
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(a) Total dose. 

Figure 5. Doses in tissue behind 20 g/cm 2 of iron shield to normal incidence of a Webber 
solar-flare proton spectrum. 100-MV rigidity. 
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(b) Total dose equivalent. 
Figure 5. Continued. 
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(c) Primary dose. 
Figure 5. Continued. 
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(e) Total secondary and heavy-ion-recoil dose. 
Figure 5. Continued. 
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x, g/cm 2 


(a) Dose. 

Figure 6. Dose in tissue behind various thicknesses of aluminum shield to normal incidence of 
February 1956 solar-flare event. 
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